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The GARCH algorithm is the most renowned generalisation of Engle's original proposal for 
modelising returns, the ARCH process. Both cases are characterised by presenting a time depen- 
dent and correlated variance or volatility. Besides a memory parameter, 6, (present in ARCH) 
and an independent and identically distributed noise, u>, GARCH involves another parameter, c, 
such that, for c = 0, the standard ARCH process is reproduced. In this manuscript we use a 
generalised noise following a distribution characterised by an index q n , such that q n — 1 recovers 
the Gaussian distribution. Matching low statistical moments of GARCH distribution for returns 
with a q- Gaussian distribution obtained through maximising the entropy S q — — ^z"f~S basis of 
nonextensive statistical mechanics, we obtain a sole analytical connection between q and (b,c,q n ) 
which turns out to be remarkably good when compared with computational simulations. With this 
result we also derive an analytical approximation for the stationary distribution for the (squared) 
volatility. Using a generalised Kullback-Leibler relative entropy form based on S q , we also analyse 
the degree of dependence between successive returns, zt and zt+i, of GARCH(1, 1) processes. This 
degree of dependence is quantified by an entropic index, q op . Our analysis points the existence of 
a unique relation between the three entropic indexes q° v , q and q n of the problem, independent of 
the value of (b, c). 



I. INTRODUCTION 



The study of time series plays an important role in 
science due to their ubiquity in both natural and ar- 
tificial systems. They can be found, e.g., in geoseis- 
mic (earthquakes), meteorological (El Nino), physiolog- 
ical (electroencephalographic profiles) or financial phe- 
nomena [1-5]. They comprise a sequentially ordered set 
of random variables, correlated or not, following a cer- 
tain probability function. For the uncorrelated case, the 
most perceptive is to consider the time series as a succes- 
sion of values that are associated to the same probability 
distribution, like it occurs for the ordinary random walk, 
where probability of a certain jump value is constant in 
time. This kind of process is defined as homoskedastic. 
However, there are phenomena for which the probabil- 
ity distribution associated to the stochastic variable at 
some time step t depends explicitly on t and these are 
named heteroskedastic. A simple way to obtain a het- 
eroskedastic process is to consider the same probability 
functional for all times, but with a varying second-order 
moment (or width). In this sort of stochasticity we can 
include financial time series, namely return time series, 
where second-order moment time dependence is a fea- 
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ture more than well-known [6]. Aiming to mimic this 
type of systems, Engle introduced in 1982 the autoregres- 
sive conditional heteroskedasticity (ARCH(s)) (s will be 
defined later on) process [7], which is considered a land- 
mark in finance, comparable to the Black-Scholes equa- 
tion [8], because of its wide use [9,10]. Albeit its ex- 
traordinary importance, Engle's model presents, in some 
situations, implementation difficulties and above that, is 
not able to fully satisfactorily reproduce the empirically 
observed power-law like correlation decay of volatility. 
These failures lead T. Bollerslev to generalise it defining 
the GARCH(s,r) [11] (G stands for generalised). Due 
to its financial cradle these models are not well-known 
in physics, nevertheless they can be very useful to illus- 
trate many traditional physical problems (see, e.g., [12]). 
In this manuscript we give sequence to the ansatz pre- 
sented by us in a previous work [13] for the ARCH (1) 
process, generalising it to the GARCH(1, 1). Moreover, 
we detail the physical justification for why these models 
(which present time dependent variance), accommodate 
well within the current nonextensive statistical mechan- 
ical theory. We also present an analytical form for the 
distribution for the (squared) volatility and analyse the 
degree of dependence between successive returns. The 
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manuscript is organised as follows: In Sec. II we intro- 
duce the GARCH(s,r) process and some of its prop- 
erties; In Sec. Ill A we present our connection between 
GARCH(1, 1) and nonextensive entropy. In addition, 
along the lines of superstatistics [15], we derive the dis- 
tribution for the (squared) volatility in Sec. IIIB. In 
Sec. IV, applying a generalised Kullback-Leibler relative 
entropy we analyse the degree of dependence between z t 
and Zt+i elements of a GARCH(1, 1) time series and its 
relation with non-Gaussianity. Some final comments are 
made in Sec. V. 



II. THE GARCH MODEL 

As settled by Engle [7] , we will define an autoregressive 
conditional heteroskedastic time series as a discrete time 
stochastic process, zt, 

z t = a t u t , (1) 

where ui t is an independent and identically distributed 
random variable with null mean and unitary variance, 
i.e., (tat) = and (wf) = 1. The quantity at, named 
volatility is time varying, positive defined and dependent 
of the past values of the return z t . According to its def- 
inition, the process presents mean zero, is uncorrelated 
((z t z t >) ~ 5 t t') and has a conditional variance, of, that 
evolves in time. 

In his original paper Engle [7] suggests a possible ex- 
pression for <7j defining it as a linear function of past 
squared values of z t known as ARCH(s) linear process, 




FIG. 1. Typical GARCH (1, 1) time series obtained for a 
Gaussian noise, q n = 1. In (a) we present time series which 
correspond to pure Gaussian (□) and ARCH (1) (o). In (b) 
the introduction of parameter c (o) increases the probability 
for larger values of \z t \, thus increasing tails in P (z). 



{a,bi ^ 0). 



(2) 



With the aim of solving the weak points of the ARCH 
process that were referred in Sec. I, it was introduced 
the linear GARCH(s,r), which presents a more flexible 
structure for the functional form of a\ decreasing the di- 
rect influence of z on of (for details see [11]): 



+ J2 biZ t-i + J2' 



(a,b i ,c l ^ 0) , (3) 



i=i 



Like its predecessor, the GARCH(s, r) model also cap- 
tures the recognised tendency for volatility clustering (ev- 
ident in financial time series) and is very similar to inter- 
mittent fluctuations in turbulent flows [14]: large(small) 
values of z t are usually followed by large(small) values. 
However, due to the arbitrary sign of oj t , it can be ver- 
ified that, although (z t z t >} ~ 5tt' , {\zt\ \zr\) is not pro- 
portional to 5 t t>- For a — (Vi), GARCH(s,r) straight- 
forwardly reduces to the linear ARCH(s). See Figs. 1, 
2 and 3 for typical realisations. 
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FIG. 2. Time dependence of volatility a for the GARCH 
processes presented in Fig. 1. Here is clearly visible the dif- 
ference between c = (ARCH) and c / (GARCH). For 
the same values of b we are able to obtain greater values of a, 
thus leading to fatter tails in P (z). 
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FIG. 3. In (a) correlation between returns z(t + T) and z(t) 
vs. time lag T . All values, except for T = 0, are at noise 
level. In (b) correlation between absolute values of returns 
which present a decay of time lag T. The run analysed was 
the same presented in Fig. 1 

Let us focus on the simplest and most used of the 
GARCH processes, namely the GARCHil, 1) one. In 
this case, 



2 ■ i 2 i 2 

<r t =a + bz t _ 1 +ca t _ 1 , 



(4) 



so that the process is completely defined when a, 6, c, 
and the noise nature, P n (u>t), are specified. 

Combining Eqs. (1) and (4) we obtain the n th order 
moment for the stationary P (z) distribution, particu- 
larly the second, 



- 2 - <^ 2 > = <- 2 > = r 

and the fourth moment, 



(b + c<l), (5) 



/ 4\ 2/4\ 1 -|- 6 -|- C 

\ z */= a \ u */ (1-6-c) (l-26c-c 2 -6 2 (uf))- 



(6) 



Let us now assume, for simplicity and without lack of 
generality, a GARCH (1,1) process that generates time 
series with unitary variance, i.e., a 2 — 1, which imposes 
a = 1 — b — c. Now, for this process, the fourth moment 

(x 4 ) 

is numerically equal to the kurtosis k x = ' ' a , and thus 
we get, 

(4) 



k z = ku 1 + y 



i 



1-c 2 -2bc-b 2 k u 



(7) 



where c 2 + 2 b c + b 2 k u < 1. It is clear from Eq. (7) that 
GARCH(l, 1) generates distributions P(z) with fatter 



tails than those of the noise u>t, and fatter also than the 
ones obtained with ARCH(\) with the same b [16]. It 
is also possible to see that parameter c is only useful 
for b 7^ 0, otherwise GARCH(1, 1) process reduces to 
constructing stationary probability distributions with a 
kurtosis k z equal to k w . 



III. THE ANSATZ CONNECTING GARCH AND 
NONEXTENSIVE STATISTICAL MECHANICS 

A. Stationary distribution for returns 

We shall now establish a connection — physically mo- 
tivated in the next subsection — between the present 
stochastic process and the current nonextensive statis- 
tical mechanical framework, based on the entropic form 
[17], 



I - 1^ [p(z)] q dz 



(8) 



This entropy is currently referred to as nonextensive be- 
cause it is so for independent subsystems. It can however 
be extensive in the presence of scale-invariant correla- 
tions [18-21]. The associated statistics has been suc- 
cessfully applied to phenomena presenting some kind of 
scale-invariant geometry, like in low-dimensional dissi- 
pativc and conservative maps [22], anomalous (corre- 
lated) diffusion [23], turbulent flows [24], Langevin dy- 
namics with fluctuating temperature [15,25,26], long- 
range many-body classical Hamiltonians [27], among 
many others [28]. Entropy (8) constitutes a generalisa- 
tion of the Boltzmann-Gibbs (BG) one, namely Sbg — 
— J-^ p{ z ) \np(z)dz. Indeed, this celebrated expression 
is recovered as the <7 — > 1 limit of entropy (8). 

By applying the standard variational principle on en- 
tropy (8) with the constraints J^ 00 p(z) dz = 1 and 

z 2 W)\ q d*/ L + r WW ^ = ^ [29-31] (a 2 is de- 
fined as the generalised second-order moment) we obtain 



p{z) 



A 



[1 + B (q-1) z 2 }-- 



where 



B = 



a 2 (5-3g)' 



1 < 3 ) , (9) 



(10) 



/ + OC 
p( 
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z) dz, 



and 



A = 



r 


i 






' 3- 
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(11) 



(12) 
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Standard and generalised second-order moments, a 2 
and a 2 are related by a 2 — a 2 _ 3 ~ <? . In the limit 
q — * 1, distribution (9) becomes a Gaussian. If <? 



f±^(m = 1,2,3. 

l+m V 



distribution (9) recovers the Stu- 
dent's t-distribution with m degrees of freedom; if q = 
(n — 3, 4, 5...), it recovers the so-called r-distribution 
with n degrees of freedom [32]. 

Defining the q- exponential function as, 



expq(x) = [1 + (1 — q)x] l -i (expi(x) = e x ) 

we can rewrite the above distribution as follows, 

p(z)=Ae- B *\ 

from now on referred to as q-Gaussian. The fourth mo- 
ment of p(z) is, 

_ 2 ^2 3q- 5 



(13) 



(14) 



(15) 
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Let us now propose the ansatz p{z) ~ P(z) with a 2 = 1. 
Specifically, we will impose the matching of Eqs. (7) and 
(15). We assume that the noise u t follows the generalised 
distribution 



Pn(w) 



1 



5-3q„ 



1 2 



(16) 



defined by the index q n ; its variance equals unity, and 
A Qn is uniquely determined through normalization. We 
are then able to establish a relation between parameters 
b, c and indices q n and q: 



y / (q-qr l )[(5-3q n )(2-q)-c 2 (5-3q)(2-q~] 
(5-3?„)(2-g) 

c(g-gn) 

(5-3g„)(2- 9 ) • 



(17) 



For 6 = c = we have q = q n . Naturally, the con- 
nection indicated in Eq. (17) and illustrated in Fig. 4 
for typical values of q n reduces, for c = 0, to the one 
obtained [13] for the ARCH (1) model, namely b = 
V(<l-<ln)/{(5-3q n )(2-q)}, hence q = [q n + 26 2 (5 - 
3q n )]/[l + b 2 (5-3q n )}. 

To verify the above ansatz we generated, for typi- 
cal values of (6, c, q n ) and using an algorithm based on 
Eqs. (1) and (4), a set of GARCH time series. Then 
we computed the corresponding probability density func- 
tions and compared them with the histograms (with 
any adequately chosen interval 5) associated with the q- 
Gaussian distribution with q satisfying the ansatz. We 
compared then the numerical probability density func- 
tions (PDFs) with H{z) 

h{z) = f:::;* P (x) dx = ^^^^ 



1-1) (5-2 zY 
L (3 qr-5) 



+ (d + Zz) 2 i<l ^2'— >2> 4(3q-5) )) 



(18) 



where 2-Fl is the hypergeometric function. 
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FIG. 4. Diagram (q,b,c) for q n — 1 in a) and q n — 1.2 
in b). In b) we can see that the greatest allowed value of b 
(c = 0) is b = -L= ~ 0.488. 



As can be seen in Figs. 5 and 6, the accordance be- 
tween stationary PDFs and the PDFs obtained by using, 
in Eq. (9), the value of q satisfying Eq. (17), is quite sat- 
isfactory. In the captions of Figs. 5 and 6, we present the 
values of the \ 2 error function, 
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Another way to evaluate the slight discrepancy between 
P (z) and p(z) (or H(z)) is to compute the percentual 
error in the sixth-order moment between PDFs. The 
results presented in Tab. I show that discrepancies are 
never larger than 3%, which in practice can be consid- 
ered negligible. It is interesting to notice the remarkable 
agreement at least down to p (z) = 1CT 6 (typical limit 
used in finance, see Ref. [4]). 




FIG. 5. PDFs for a q n = 1 noise and typical values ol (6; c) 
pair, a) (0.1; 0.1), q = 1.021 ( X 2 = 2.35xl0~ 9 ); b) (0.1; 0.88), 
q = 1.287 (x 2 = 4.59 x 10~ 10 ); c) (0.4; 0.1), q = 1.26 
(X 2 = 2.44 x 10~ 9 ); d) (0.4; 0.4), q = 1.38 (x 2 = 3.22 x 10~ 7 ). 




FIG. 6. PDFs for a q n = 1.2 noise and typical values ol 
(6; c) pair, a) (0.1; 0.1), q = 1.211 (x 2 = 8.67 x 10~ 10 ); 
b) (0.1; 0.5), q = 1.221 (x 2 = 6.01 x 10~ 10 ); c) (0.3; 0.25), 
q = 1.310 (x 2 = 8.11 x 10~ 9 ); d) (0.3; 0.45), q = 1.35 
(X 2 = 7.36 x 10~ 9 ). 



B. Stationary distribution for (squared) volatility 

The good results of the proposal presented in the pre- 
vious sub-section have their explanation in the main fea- 
ture of the ARCH processes, namely the temporal depen- 
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dence of the variance. This fact allow us to think of them 
along the same lines of Wilk and Wlodarczyk [25] and of 
Beck [26] . They essentially considered a system following 
a Langevin equation with fluctuating temperature, gene- 
sis of the superstatistics (statistics of statistics) recently 
advanced by Beck and Cohen [15]. This approach was 
developed to treat non-equilibrium systems composed, 
for instance, of smaller cells obeying BG statistics with a 
distribution Pbg °^ e~@ E , E being the energy. The long- 
term stationary state that presents a spatio/temporaly 
fluctuating temperature following a distribution f((3). In 
the long-term, the probability density for the nonequi- 
librium system comes from BG statistics associated with 
the small cells that are averaged over the various (3, i.e., 



stationary 



(E) = f /(/?) P BG (E)dp 



(19) 



Instead of defining the intensive fluctuating parameter, 
P, as the inverse temperature, we will define it as an in- 
verse second-order moment, [3 a = j^. For long times, 
we can study GARCH as a stationary diffusion process, 

composed by t increments, each of them of size z t . The 

t 

total value is given by x t = X) z i' wnere each increment 

i=i 

follows a certain q n — Gaussian, with a second-order mo- 
ment, cr 2 , that is associated with a probability distribu- 
tion, p a (a 2 ). The stationary probability distribution p(z) 
is thus given by 



p{z) 



-f 

Jo 



p a {a 2 ) p(z\a 2 ) d(a 2 ), 



(20) 



where p(z\a 2 ) is the conditional probability of having a 
value z for the return given a value a 2 . The homoskcdas- 
tic case corresponds to p a (cr 2 ) = S (a 2 — cr 2 ). Let us 
focus, for now, on the case q n = 1 (Gaussian noise) with 



p(z\a 2 ) = 



V2ira'- 



. e 2 " i 



(21) 



In their proposal Beck and Cohen [15] showed that, if the 
intensive parameter /3 is associated with a Gamma dis- 
tribution, then the macroscopic non-equilibrium steady 
state follows exactly a q— exponential distribution (see 
also [25,26]). Following the reverse line and assuming 



p{z) = P ( 



stationary 



~ P(z), we are lead to a Gamma dis- 



tribution in f3 a . In other words, 



Pa \o 



ex P(-7^) ( g2 ) 
(2 n) X T[A] 



-l-A 



where 



and 



1 1 



q- 1 2 

1-g 
a 2 (3g-5) ' 



(22) 



(23) 



(24) 



As can be seen in Fig. 7, the ansatz gives also a quite 
satisfactory description for the probability distribution 
in the (squared) volatility, corroborating the connection 
between the ARCH class of processes, nonextensive sta- 
tistical mechanics, and superstatistics. The analytic ex- 
pressions (22), (23) and (24) can be very useful in ap- 
plications related, among others, to option prices [4,33] 
where volatility forecasting plays a particularly impor- 
tant role [34,35]. 
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FIG. 7. The symbols in black represent the cumulative 
probability distribution, (a 2 ) obtained from numerical 
simulation for a Gaussian noise with b — c = 0.4 . The gray 
line represents the same distribution as given by Eq. (22) with 
(k, A, a 2 ) = (0.444,2.125, 1) satisfying Eqs. (23) and (24). 

For the case of a (7 n -Gaussian noise (with q n > 1) simi- 
lar arguments can be used. However, the achievement of 
an analytical solution for Eq. (20) which generalises the 
approach above is not a trivial task. 

For the conditional probability p(z \a 2 ) we know a sat- 
isfactory answer [36] , namely 



p(z\a 2 ) 



A 



[l + B ( , B ,„ a) (q n -l) z 2 }^ 



(25) 



But the same does not happen with p a (er 2 ) . Assuming 
a continuous approach in q n , a good ansatz for describing 
Pa (er 2 ) is 



p a ((T 2 ) OC (cr 2 ) 1 A CXp^ 

=(- 2 )- 1 - A (l+ 



1 



2 K G 2 

q<j - 1 1 \ " " 

2 k cr 5 



(26) 



where q„ is an index which depends on q n , i.e., q a {q n ) 
such that q a (1) = 1. For large cr 2 , 



Pa (a 2 ) - (a 2 ) 



-l-A 



(27) 



The integral J °° a 2 p a {a 2 ) d(a 2 ) should equal the mean 
value a 2 = 1 _ a b _ c - This yields, 
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1 + A - A q a 

2 k (A - 1) 



= a 



and 



/oo />oo 
z 4 / p a (a 2 ) p (z | a 2 ) d (a 2 ) dz = 3 (a 2 ) 
-oo JO 



(28) 



2 3 g — 5 
5(7-7' 

(29) 



From Eq. (27) and by adjusting numerically the curves 
for the cumulated probability distributions, (er 2 ) , we 
were able to determine A, and then q a and k from 
Eqs. (28) and (29). The procedure appears to be valid for 
values of q n close to unity. From the analysis of some val- 
ues of q n we verified that the value of q a equals 1 for every 
q n considered. Figures 7 and 8 confirm that our pro- 
posal produces a satisfactory approximation when com- 
pared with numerical simulations, particularly for large 
values of volatility (which are, in turn, responsible for 
the large returns). For the q n — Gaussian noise case, al- 
though some discrepancy exits for small er 2 , the tail is 
remarkably good. 




FIG. 8. The line in black represent the cumulative proba- 
bility distribution, P^ (a ) obtained from numerical simula- 
tion for a q n — Gaussian noise (q n — 1.15) with (b, c) = (0.5, 0). 
The gray line represents the same distribution as given by Eq. 
(26) with (k, A, a 2 ) = (0.365, 2.371, 1), and q„ = 1. 



IV. DEGREE OF DEPENDENCE BETWEEN 
SUCCESSIVE RETURNS 

As stated in Sec. II, stochastic variables, {z t } 1 in a 
G ARC 'H process are uncorrelated (see Eq. (1)). How- 
ever, if we combine Eqs. (1) and (3), we can verify 
that they are not independent. More specifically, for 
GARCH(1, 1) (a 2 = 1), we have 



zt 



'l + b(z 



\ w t-i 



(30) 



One of the possible measures of the dependence between 
the (z t ; Zt—i) stochastic variables consists in using a gen- 
eralised form of the Kullback-Leibler relative entropy 
[37], namely 



I q >(p 1 ,p 2 ) = - J P\{u) Lry 



Pi( u ). 



du , 



(31) 



where ln g / (a;) 



- Yz^p 1 ((/-logarithm). In the limit 
q' — » 1, I q i recovers the standard Kullback-Leibler form 
[38] . This generalised relative entropy equals zero when- 
ever p2(u) = pi(u), and has the same sign as q' other- 
wise. With these properties we can use I q i (with q' > 0) 
as a way to compute the "distance" , in probability space, 
from p2 (u) to pi (u) . Assume that u is a two dimensional 
random variable u = (x, y) so that pi(x, y) represents the 
joint distribution of (x,y), and p 2 (x,y) = /ii(x) /i 2 (y), 
where h\(x) = J pi(x,y) dy and ti2{y) = J pi(x,y) dx 
are the marginal distributions. Random variables x and 
y are independent if pi(x,y) = P2(x,y). With the func- 
tional I q i we can measure the degree of dependence via 
distance between probabilities pi(x, y) and p 2 (x, y). 




FIG. 9. Dependence criterion R q i vs. q' for various 
GARCH(1, 1) process with typical (q n ,b,c) triplets. In 
(a) q n — 1 and (6, c) have been chosen as follows: l-(0, 0), 
2-(0.05,0),3-(0.1,0.2),4-(0.15,0),5-(0.2,0),6-(0.2,0.2),7-(0.25,0), 
8-(0.3,0), 9-(0.4,0), 10-(0.4,0.1), ll-(0.4, 0.2), 12-(0.4, 0.4), 
13-(0.5,0), 14-(0.2, 0.688), 15-(0.35,0), 16-(0.1,0). In (b) 
q n = 1.2 and (b, c) have been chosen as follows: l-(0, 0), 
2-(0.1,0), 3-(0.1,0.1), 4-(0.15,0), 5-(0.2,0), 6-(0.25,0), 
7-(0.3,0.1), 8-(0.377,0), 9-(0.3, 0.45), 10-(0.48, 0.0). The in- 
sets contain the derivative dR q i/dq' (numerically obtained) 
for the first four curves as mere illustration. 
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For this case, u = (x, y) and q' > 0, it was shown [40] 
that, besides a lower bound I q i = (total independence 
of x and y) , I q > presents, for every value of q , an upper 
bound (complete dependence between x and y) given by, 

i q r AX (p 1 ,p 2 ) = - // \pi(x,y)} q ' ^(x) 

+ (1 - q')[\n ql hr{x))[\n q , h 2 (y)}} dx dy . (32) 



Dividing Eq. (31) by (32) we define 



Rqi = 



jMAX 



G [0,1] • 



(33) 



This ratio can be used as a criterion for measuring the 
degree of dependence between random variables. Indeed, 
it presents an optimal q', noted q op [39], for which the 
sensitivity of R q i is maximal (q op corresponds to the in- 
flexion point of R q '(q')). Higher (lower) values of q op 
represents lower (higher) degree of dependence [40] . Tak- 
ing x = z t and y — z t -\ generated from Eq. (30), and 
applying Eq. (33), we obtained the curves presented in 
Fig. 9 for typical values of (b, c, q n ). For each set we 
determined q op and plotted it versus q obtained from 
Eq. (17): see Fig. 10. For both noises that have been 
illustrated, q op monotonically decreases with q. For fixed 
q n , this (decreasing) curve does not depend on the values 
of (&, c) . An illustration of this independence is indicated 
in Fig. 10 by using different pairs (b, c) that give the same 
q. This fact suggests the existence of a relation between 
non-Gaussianity (represented by q), the degree of depen- 
dence quantified with q op , and the noise index q n . This 
triangular relation {q,q op ,q n ) is analogous to another one 
which could relate the sensitivity, relaxation and station- 
arity in weakly chaotic systems such as those in which 
long-range interactions are assumed [41]. 



2.50^ 



2.25- 



2.00- 



1.75 



1.50 J 



D 1n = 1 

o q°=1.2 

1 n 



b = 0.4, c = 
b = 0.2, c = 0.6 
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q 

FIG. 10. Plot of q op vs. q for typical (q n ,b,c) triplets. The 
arrow points two examples which were obtained from differ- 
ent triplets, and nevertheless coincide in what concerns the 
resulting point (q,q op ). 



cess, which is characterised by the fact that it exhibits 
a time dependent (and correlated) volatility. This at- 
tribute makes this system similar to those exhibiting fluc- 
tuations in some intensive parameter (e.g. temperature), 
which is at the basis of superstatistics. Inspired by the 
intimate connection between superstatistics and nonex- 
tensive statistical mechanics, we obtained an expression 
linking the dynamical parameters (b, c) and the noise na- 
ture, q n , with the entropic index q characterising the 
stationary distribution for the associated GARCH(1, 1) 
process. After numerically testing the validity of the 
approach for the return distributions, we derived ana- 
lytical expressions for the squared volatility stationary 
distribution in the presence of g„-Gaussian noise. The 
results are satisfactory for q n ~ 1. Then, using the q- 
generalised Kullback-Leibler relative entropy, we quanti- 
fied the degree of dependence between successive returns. 
This analysis led to an entropic index q op which is optimal 
in the sense that the ratio (33) exhibits maximal sensitiv- 
ity We then verified the existence of a direct relation be- 
tween q op , the non-Gaussianity, q, and the nature of the 
noise, q n . An interesting property emerged, namely that, 
whatever be the pair (b, c) that results in a certain q for 
the stationary distribution, one obtains the same value of 
q op . Consequently, the time series will present the same 
degree of dependence. The connection between various 
entropic indices reminds us of the dynamical scenario 
within which nonextensive statistical mechanics is for- 
mulated. Indeed, various entropic indices emerge therein, 
which coalesce onto the same value q = 1 if ergodicity is 
satisfied. In the present context, this connection can be 
analysed as follows. Due to the dynamical characteristics 
of GARCH(1, 1), non-Gaussian distribution for returns 
(q 7^ 1) comes from temporal dependence on its second- 
order moment which is a self-correlated variable. This 
correlation, sign of memory in the process, is the respon- 
sible for the breakdown of independence [19] between z t 
and Zt+i, which in turn reflects on the entropic index 
q op . The connection between q and q op opens the door 
to the establishment of a relation between the topologies 
of phase space and probability space. A careful analy- 
sis of other kind of systems (e.g. long-range hamiltonian 
models, Langevin-like dynamics obeying generalised cor- 
related Fokker-Plank equations) should give a deep in- 
sight onto this question. 
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V. CONCLUDING REMARKS 



In this article, we presented a study of the stationary 
statistical features of the GARCH(l, 1) stochastic pro- 
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q n 


b 


c 


\ / n.n.m.p.ri.r.a.1 


\ I am.xa.t7. 


error (%) 


1 


0.1 


0.2 


13.97 


14.05 


0.60 


1 


0.1 


0.88 


66.75 


65.46 


1.93 


1 


0.4 


0.1 


61.45 


62.33 


1.44 


1 


0.4 


0.4 


591.91 


517.67 


1.75 


1.2 


0.1 


0.1 


49.41 


49.08 


0.67 


1.2 


0.1 


0.5 


55.93 


55.43 


0.89 


1.2 


0.3 


0.25 


181.48 


179.44 


1.13 


1.2 


0.3 


0.45 


1416.41 


1455.37 


2.75 



TABLE I. Percentual error in the sixth-order moment be- 
tween numerical and ansatz PDFs presented in Figs. 5 and 
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